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I review recent progress in the calculation of inclusive quarkonium decays in the framework of QCD nonrela- 
tivistic effective field theories and in relation to the experimental measurements. 



1. INTRODUCTION 

Inclusive quarkonium decays, more precisely 
inclusive annihilation decay rates of heavy 
quarkonium states into light hadrons (hadronic 
decays) and photons and lepton pairs (electro- 
magnetic decays), were among the earliest calcu- 
lations of perturbative QCD. There, it was as- 
sumed that the decay rate of the quarkonium 
state factored into a short distance part, calcu- 
lated as the annihilation rate of the heavy quark 
and antiquark and given in terms of as{m), and 
a long distance nonperturbative part given in 
terms of the quarkonium wave function (or its 
derivatives) evaluated at the origin. Explicit cal- 
culations at next to leading order in as in per- 
turbation theory for S- and P-wave decays sup- 
ported the factorization assumption which could 
not however be proved on general grounds for 
higher orders of perturbation theory. Indeed, in 
the case of P-wave decays into light hadrons, it 
turned out that at order the factorization was 
spoiled by logarithmic infrared divergences [1]. 
The same problem appeared in relativistic correc- 
tions to the annihilation decays of S'-wave states 

[!]• 

This problem has been solved [2,3] and new pre- 
dictions have been obtained [4,5] with the intro- 
duction of non-relativistic effective field theories 
(EFTs) of QCD, that has put our description of 
these systems on the solid ground of QCD. 

The reason for which the EFT approach is so 
successfull for heavy quarkonium is the follow- 
ing. Heavy quarkonium, being a non-relativistic 
bound state, is characterized by a hierarchy of 
energy scales to, mv and tou^, where to is the 



heavy-quark mass and v <^ 1 the relative heavy- 
quark velocity. A hierarchy of EFTs may be con- 
structed by systematically integrating out modes 
associated to these energy scales in a matching 
procedure that enforces the complete equivalence 
between QCD and the EFT at a given order of 
the expansion in v and ag. 

Integrating out degrees of freedom of energy 
m, which for heavy quarks can be done perturba- 
tively, leads to non-relativistic QCD (NRQCD)[2, 
3] . This EFT still contains the lower energy scales 
as dynamical degrees of freedom. In the last 
years, the problem of integrating out the remain- 
ing dynamical scales of NRQCD has been ad- 
dressed by several groups and has reached a good 
level of understanding (a list of references can be 
found in [6]). The EFT obtained by subsequent 
matchings from QCD, where only the lightest de- 
grees of freedom of energy mv^ are left dynami- 
cal, is called potential NRQCD, pNRQCD [4,5]. 
This EFT is close to a quantum-mechanical de- 
scription of the bound system and, therefore, as 
simple. It has been systematically explored in the 
dynamical regime Aqcd ^ ^tiw^ in [5,7,8] and in 
the regime mv'^ <C Aqcd ^ "mv in [5,9,10]. The 
quantity Aqcd stands for the generic scale of non- 
perturbative physics. 

The EFT approach made it possible, in the case 
of several observables, among which the inclusive 
decay widths, to achieve a rigorous factorization 
between the high-energy dynamics encoded into 
matching coefficients calculable in perturbation 
theory and the non-pcrturbative QCD dynamics 
encoded into few well-defined matrix elements. 
Systematic improvements are possible, either by 



calculating higher-order corrections in the cou- 
pling constant or by adding higher-order opera- 
tors. In this way the prediction of inclusive decays 
are put in direct relation to QCD and the theoret- 
ical uncertainties may be consistently estimated. 

Prom the experimental side new data have re- 
cently been produced for heavy- quarkonium ob- 
servables. Measurements relevant to the de- 
termination of heavy-quarkonium inclusive de- 
cay widths have come from Fermilab (E835)[ll], 
BES[12], CLEO[13,14] and Belle[15] ^ and de- 
mand accurate QCD theoretical predictions. On 
the other hand the inclusive decays of heavy 
quarkonium may provide competitive information 
on as at the scale m once the theoretical and ex- 
perimental errors are under control. 

In the following I will review recent progress 
in our theoretical understanding of inclusive and 
electromagnetic heavy-quarkonium decays. I will 
recall the NRQCD factorization results in Sec. 2 
and the further simplification achieved within the 
pNRQCD factorization in Sec. 3. The presented 
pNRQCD formulas apply to quarkonia that fulfil 
mv > Aqcd » rnv'^. 



2. NRQCD 

NRQCD is the EFT obtained by integrating 
out the hard scale m. This m being larger than 
the scale of non-perturbative physics, Aqcd, the 
matching to NRQCD can be done order by order 
in as- Hence, the NRQCD Lagrangian can be 
written as a sum of terms like fnOn"^ /m'^"~^, 
ordered in powers of as and v. More specifically, 
the Wilson coefficients are series in as(m) and 
encode the ultraviolet physics that has been in- 
tegrated out from QCD. The operators oi'^"'' of 
dimension dn describe the low-energy dynamics 
and are counted in powers of v. Heavy quarko- 
nium inclusive decays arc controlled by the imag- 
inary part of the NRQCD Hamiltonian, i.e. the 
imaginary part of the Wilson coefficients of the 

4-fermion operators {On"'^ = i)^ KnXX^ ^'ni') 
the NRQCD Lagrangian. The NRQCD factoriza- 
tion formula for quarkonium (H) inclusive decay 



widths into light hadrons (LH) reads [3] 
r(H - LH) = 5] (H|Vti^„XX^<V'|H).(l) 



The 4-fermion operators are classified with re- 
spect to their rotational and spin symmetry (e.g. 
0(2^+1 Sj), 0(2^+iPj), ...) and of their colour 
content (octet, Os, and singlet, 0\, operators). 
Singlet operator expectation values may be eas- 
ily related to the square of the quarkonium wave 
functions (or derivatives of it) at the origin. 
These are unknown non-perturbative parameters. 

Let us make a concrete example by considering 
the P-wave inclusive decays. 

In NRQCD the P-wave inclusive decay width 
for the S = [h) and S = \ (x) quarkonimn 
states is given at leading (non- vanishing) order in 
V (which is mv^) by [3]: 



T{h LH) 



9Im/i(iPi) 



Pp 



(2) 



r(xj ^ LH) 



9Im/i(3Pj) 



7rm* 



Pp 



+ 



^^"'^f^''' {h\OsCSo)\h)MJ = 0, 1, 2 (3) 



'^For updates see also [25]. 



where R'p is the derivative of the radial P-wave 
function at the origin. We stress that, accord- 
ing to the power counting of NRQCD, the octet 
contribution {h\OsCSo)\h) is as relevant as the 
singlet contribution [3]. This octet contribution 
reabsorbs the dependence on the infrared cut-off 
/X of the Wilson coefficients Im /i (P) solving the 
problem mentioned in the introduction. From the 
above equations, we see that in NRQCD the 8 
P-wave bottomonium states (IP, 2P) and the 4 
P-wave charmonium states (IP), which lie under 
threshold, depend at leading order in the veloc- 
ity expansion on 6 non-perturbative parameters 
(3 wave functions + 3 octet matrix elements). 
The inclusive decays of the heavy quarkonium 
(either hadronic or electromagnetic) are usually 
considered up to, and including, NRQCD matrix 
elements of 4-fermion operators of dimension 8. 
This means to consider the ©(l/m^, l/m^) local 
4-fermion operators of the NRQCD Lagrangian 



and the decay rate up to order mv^ in the v count- 
ing. If we consider that in the bottomonium sys- 
tem in principle 14 5*- and P-wave states lie below 
threshold {T{nS) and r]b{nS) with n = 1,2,3; 
hb{nP) and Xb.j{nP) with n = 1,2 and J = 
0, 1,2) and that in the charmonium system this 
is the case for 8 states {■tp{nS) and rjdnS) with 
n = 1,2; K{IP) and Xcj{^P) with J = 0,1,2), 
all the bottomonium and charmonium S- and P- 
wave decays into light hadrons and into photons 
or e+e~ are then described in NRQCD up to 
by 46 unknown NRQCD matrix elements (40 for 
the 5- wave decays and 6 for the P-wave decays) , 
where we have already used spin symmetry and 
vacuum saturation [3,10]. These matrix elements 
have to be fixed either by lattice simulations [16] 
or by fitting the data [17]. Only in the spe- 
cific case of matrix elements of singlet operators 
does NRQCD allow an interpretation in terms of 
quarkonium wave functions and one can resort to 
potential models. 

It has been discussed, in particular in [18] (but 
see also [19]), that higher-order operators, not 
considered in the above formulas, can be numeri- 
cally quite relevant. This may be the case partic- 
ularly for charmonium, where ~ 0.3, so that 
relativistic corrections arc large, and for P-wave 
decays where the above formulas provide, indeed, 
only the leading-order contribution in the velocity 
expansion. Besides this, precise theoretical pre- 
dictions are also hampered by uncertainties in the 
NRQCD matrix elements and large corrections in 
NLO in Us- The convergence of the perturbative 
series of the four-fermion matching coefficients is 
indeed often bad (for examples see [22]). A so- 
lution may be provided by the resummation of 
the large contributions in the perturbative series 
coming from bubble-chain diagrams. This analy- 
sis has been successfully carried out in some spe- 
cific cases in [20] but a general treatment is still 
missing, in particular in the case of P-wave de- 
cays. 

Here we have counted NRQCD matrix ele- 
ments by their dimensionality only. In fact the 
power counting of the NRQCD matrix elements 
is an open issue. In the standard NRQCD power 
counting, the octet matrix elements are 0{v^) 
suppressed for S-wave decays if compared with 



the leading order. This is not so within pNRQCD, 
assuming the counting Aqcd ^ mv, they would 
only be C'(?;^)-suppressed [10]. This is poten- 
tially relevant for T{V LH) since lxafi{^Si) 
is C'(as('rri))-suppressed with respect to Im/8(S') 
[10]. In other words, the octet matrix element ef- 
fects could potentially be much more important 
than usually thought for these decays. 

In the next section we will show that in the 
framework of pNRQCD it is possible to achieve a 
noticeable reduction in the number of the nonper- 
turbative parameters and thus to formulate new 
predictions with respect to NRQCD. 

3. pNRQCD 

Pushing further the EFT programme for non- 
relativistic bound states, further simplifications 
occur if we integrate out also soft degrees of free- 
dom. pNRQCD is the resulting EFT. We will 
consider pNRQCD under the condition Aqcd ^ 
mv^ . Then, two situations are possible. First, the 
situation when mv » Aqcd ^ mv"^ . In this case 
the soft scale mv can be integrated out pertur- 
batively. This leads to an intermediate EFT that 
contains singlet and octet quarkonium fields and 
ultrasoft gluons as dynamical degrees of freedom. 
The octet quarkonium field and the ultrasoft glu- 
ons are eventually integrated out by the (non- 
perturbative) matching to pNRQCD [5] . Second, 
the situation when Aqcd ^ "mv. In this case the 
(nonperturbative) matching to pNRQCD has to 
be done in one single step [9]. Under the circum- 
stances that other degrees of freedom develop a 
mass gap of order Aqcd the quarkonium singlet 
field S remains as the only dynamical degree of 
freedom in the pNRQCD Lagrangian, which reads 

[5,9,10] £pNRQCD = Tr{st(«ao-W)s}, U be- 
ing the pNRQCD Hamiltonian, to be determined 
by matching pNRQCD to NRQCD. The inclu- 
sive quarkonium decay width into light hadrons 
is given by 

r(H ^ LH) = -2 Im (n, L, S, J\n\n, L, S, J), (4) 

where \n,L,S,J) is an eigenstate of H with the 
quantum numbers of the quarkonium state H. 



3.1. Matching in a 1/m expansion 

We consider first the case in which the match- 
ing between NRQCD and pNRQCD is made 
within a 1/m expansion. In this case from the 
matching we obtain schematically: 

ImH = ^^(r) Y: + {S'ir), A} (5) 

n 

X E + v^^V) v^' E iii^c^' + • • • ' 

n n 

where the imaginary part of fn are inherited 
from the 4-heavy-fermion NRQCD matching co- 
efficients, and Arn Bn, ... arc nonpcrturbativc 
operators, which are universal in the sense that 
they do not depend either on the heavy-quark 
flavour or on the spc'dfic; quantum numbers of 
the considered heavy-quarkonium state. Insert- 
ing Eq. (5) into (4) and comparing with Eq. (1), 
we see [10] that all NRQCD matrix elements, in- 
cluding the octet ones, can be expressed through 
pNRQCD as products of universal nonperturba- 
tive factors by the squares of the quarkonium 
wave functions (or derivatives of it) at the ori- 
gin. In [10] the inclusive decay widths into light 
hadrons, photons and Icpton pairs of all S'-wavc 
and P-wave states (under threshold) have been 
calculated up to 0{mv^ x (AQ^p/m^, E/m)) and 
0{mv^). A large reduction in the number of 
unknown nonperturbative parameters is achieved 
and, therefore, new model-independent QCD pre- 
dictions may be obtained. The universal non- 
perturbative parameters are all expressed only in 
terms of gluonic field-strength correlators, which 
may be fixed by experimental data or by lat- 
tice simulations. Thus at the same level of ac- 
curacy discussed before in NRQCD, S- and P- 
wavc bottomonium and charmonium decays are 
described in pNRQCD, under the dynamical as- 
sumption Aqcd ^ fnv'^ and within a. 1/m expan- 
sion matching, by only 19 nonperturbative pa- 
rameters. These are the 13 wave functions and 6 
universal nonperturbative parameters. 

This same approach may be useful also for 
quarkonium production. 

3.2. Contributions from the scale y^mAgcD 
Once the methodology to compute the poten- 
tials (real and imaginary contributions) and from 



these the inclusive decays, within a 1/m expan- 
sion in the matching has been developed, the next 
question that appears naturally is to which ex- 
tent one can compute the full potential within a 
1/m expansion in the case Aqcd ^ mv^. It has 
been recently shown [21] that new non-analytic 
terms arise due to the three-momentum scale 
y^mAqcD • These terms can be incorporated into 
local potentials (5^(r) and derivatives of it) and 
scale as half-integer powers of 1/m. Moreover, 
it is possible to factorize these effects in a model 
independent way and compute them within a sys- 
tematic expansion in some small parameters. 

These terms are due to the existence of de- 
grees of freedom, namely the quark-antiquark 
pair, with relative three-momentum of order 
lymAqcD- The on-shell energy of these degrees 
of freedom is of O(Aqcd), i.e. the same energy 
scale that is integrated out when computing the 
standard 1/m potentials, which corresponds to 
integrating out (off-shell) quark-antiquark pairs 
of three momentum of order Aqcd. Therefore 
both degrees of freedom should be integrated out 
at the same time. Note that the scale ^/m, Aqcd 
fulfils -y/mAQCD ^ Aqcd, from which it follows 
that at this scale we always are in the perturba- 
tive regime. The matching may be performed in 
two different ways depending on the two situa- 
tions mv ^ Aqcd or Aqcd < mv [21]. 

The result for the inclusive decays is the fol- 
lowing. In general, the size of the contributions 
coming from the scale ^mAQCD depends on the 
size of as (V^ Aqcd) [21]. For P dc x'.avs the 
leading effect turns out to be 0(mQ:s/-^m,AQCD) 
suppressed with respect to the leading con- 
tribution of order mv^ (we assume mag <C 

m,AQCD). For the S'-wavc decay widths the 
leading effect is O ((AQCD/m) (mag /y'm Aqcd)) 
suppressed with respect to the leading contri- 
bution of order mv^ . All the results fulfil the 
same factorization properties as those obtained 
in [10] and mentioned in Section 3.1, and thus 
the nonperturbative parameters are still encoded 
into few, only glue dependent, operators. 



3.3. Applications to P decays 

Here we show how the reduction obtained by 
pNRQCD in the numbers of unknown nonportur- 
bative factors makes new theoretical predictions 
possible in the case of P-wave inclusive quarko- 
nium decays into light hadrons [10,?]. 

In pNRQCD the P-wave inclusive decay widths 
are given at leading order in v ( which is v^) by 
[10]: 



T{h LH) = 
r(x7 ^ LH) = 



9Im^(3p,) + ^^^^£ 



(6) 
•(7) 



1 

f = -y dt t^gW{t,O)^ab{t,0;O)gE\0,O)) 

is the universal nonperturbative (only glue depen- 
dent) parameter that describes P-wave quarko- 
nium decays in pNRQCD. 

By comparing Eqs. (2) and (3) with Eqs. (6) 
and (7) we get at leading order in v the rela- 
tion between the octet matrix element of NRQCD 



Pp 



£/{18nm'^). The 



and £: {hlOsi^ So)\h) 
quarkonium-state dependence factorizes in the 
pNRQCD formulas. This allows some new pre- 
dictions with respect to NRQCD, which are syn- 
thetized by the formula (valid at LO in v): 



r(H(2S+inPj) ^ LH) 



r(H(2'S'+inPj/) ^LH) 

81Im/i(2^^+ip,7)+Im/8(2^'+iSs)f 



(8) 



81Im/i(2S'+ip^^ 



+ Im/8(2S' + lS5' 



where the left-hand side is a ratio between inclu- 
sive decay widths of P-wave quarkonia with the 
same principal quantum number n and the right- 
hand side no longer depends on n and has the 
whole flavour dependence encoded in the Wilson 
coefficients, which arc known quantities. 

In practice, the 12 P-wave quarkonium states, 
which lie under threshold, depend only, in pN- 
RQCD at leading (non-vanishing) order in the 
velocity expansion, on 4 nonperturbative param- 
eters (3 wave functions -|- 1 chromoelectric cor- 
relator £). The reduction by 2 in the number 



of unknown nonperturbative parameters with re- 
spect to NRQCD, allows us to formulate two new 
statements. In particular we can use the charmo- 
nium data to extract a determination of £, which 
in turn can be used to produce two new predic- 
tions for bottomonium, (at NLO): 



= 8.0±1.3, 



r(x6o(iP) - 
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r(x6o(2P) 
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r(x6i(2P) 
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or alternatively 






r(xw(iP) - 
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r(xw(2P) 
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r(x62(iP) - 
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r(x62(2P) 


^ LH) 
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The errors here refer only to the experimental 
errors in the charmonium inclusive decays data 
(that in turn produce an error on £) while no at- 
tempt has been done up to now to include also 
the theoretical error. Recently preliminary deter- 
minations of these two ratios have been produced 
by CLE0-III[14]: 19.3 ±9.8 for the first ratio and 
0.29 ± 0.06 for the second one. 

In Fig. 1 we display a plot of the first ratio 
of decay widths as a function of the factoriza- 
tion scale /i and at leading and next-to-leading 
order in the matching coefficients. As it should 
be, the result is stable in /j,. The obtained bands 
compare well with the first CLEO-III determina- 
tion, published after the completion of our work, 
r^^^^LH/rxM^LH = 19.3 ± 9.8. However, both 
theoretical and experimental determinations are 
still affected by large uncertainties and the large 
correction in the NLO of the matching coefficients 
should be put under control (e.g. via renormalon 
resummation) or at least considered in the errors. 

4. CONCLUSIONS AND OUTLOOK 

The progress in our understanding of non- 
relativistic efi'ective field theories makes it pos- 
sible to move beyond ad hoc phenomenological 
models and have a unified description of the dif- 
ferent heavy-quarkonium observables, so that the 
same quantities determined from a set of data 
may be used in order to describe other sets. 
Moreover, predictions based on non-relativistic 
EFTs are conceptually solid, and systematically 
improvable. In the framework of pNRQCD, for 
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figure is taken from [23]). 
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physical states that satisfy Aqcd ^ mw^, octet 
matrix elements may be expressed in terms of the 
wave function in the origin and some universal 
non-local gluon-field correlators obtaining a sig- 
nificant reduction of the nonperturbative param- 
eters. The same nonperturbative correlators en- 
ter also in the expression of the masses of some 
heavy-quarkonium states [7,24]. Difficulties still 
exist at the level of the control of higher-order 
corrections in the velocity and as expansion. In 
principle, the tools to overcome these difficulties 
already exist, so we expect relevant progress in 
the field from the coordinated effort of the heavy- 
quarkonium community in the near future [25]. 
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